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Abstract. We consider the class of long-range Hamiltonian systems first introduced 
by Anteneodo and Tsallis and called the a-XY model. This involves A classical 
rotators on a d-dimensional periodic lattice interacting all to all with an attractive 
coupling whose strength decays as r~ a , r being the distances between sites. Using a 
recent geometrical approach, we estimate for any d-dimensional lattice the scaling of 
the largest Lyapunov exponent (LLE) with A as a function of a in the large energy 
regime where rotators behave almost freely. We find that the LLE vanishes as N~ K , 
with k = 1/3 for < a/d < 1/2 and n = 2/3(1 - a/d) for 1/2 < a/d < 1. These 
analytical results present a nice agreement with numerical results obtained by Campa 
et ai, including deviations at small A. 



PACS numbers: 05.45. +b, 05.70.Fh, 02.40.-k 



Submitted to: J. Phys. A: Math. Gen. 



Letter to the Editor 



2 



It is well known that systems interacting via long-range interactions may exhibit 
pathological thermodynamical as well as dynamical behaviours. This issue has been 
very much debated recently. In particular, for systems governed by sufficiently long- 
range interactions decaying as r~ a , with a < d the Euclidean space dimension and 
r the interparticle distance, the Hamiltonian comes out to be non- extensive, that is 
the energy per particle diverges in the thermodynamic limit N —>■ oc. In order to 
study issues related to the links between non-extensivity and long-range interactions, 
Anteneodo and Tsallis introduced in reference Q a generalization of the Hamiltonian 
mean-field (HMF) ferromagnetic-like X-Y model in the form 



N „ N 



H 



2N . .^f., . r 



with d > a > 0. This is an Hamiltonian model of rotators placed at the sites of a d- 
dimensional lattice with indices i, j and is the shortest distance between them with 
periodic boundary conditions, so that the interaction between rotators % and j decays as 
the inverse of their distance to the power a. Here, the N rescaling function is introduced 
in order to get an extensive (i.e. of order N) potential V. This trick was not adopted 
in reference [[j], but then the authors had to conveniently rescale thermodynamical 
potentials by N. In the large N limit H 

Nl/d ( i u 

N = 1 + d / r d ~ x ~ a dr ~ I r (2) 

J \\n(N) iia = d 

As remarked in reference this does not entail energy additivity. In the one-dimensional 
case, Tamarit and Anteneodo M have shown numerically that the canonical caloric and 
magnetization curves for model (0) could be derived from the curves of the HMF model 
(recovered for a = 0, d = 1), and this result has been later derived analytically by 
Campa, Giansanti and Moroni |13[ . 

Recently, the stochasticity exhibited by system (JT|) has been investigated 
numerically through the computation of the largest Lyapunov exponent (LLE) as a 
function of the energy density e = H/N 0, ITSf . In the phase where particles behave 
almost freely, i.e. for e above the critical energy density e c , where the system exhibits 
a second order phase transition, the LLE has been found to scale with the number of 
rotators as N~ K where k is a so-called "universal" function of the ratio a/d with no 
dependence on the energy e nor on d. In the past, the LLE had been studied for the 
HMF model § % §[. 

For Hamiltonian models having a large number of degrees of freedom like (|1|) and 
using a geometric reformulation of the dynamics that associates trajectories to geodesies, 
Casetti, Pettini and coworkers || |H], [T1J have proposed an expression of the LLE in 
terms of the ensemble-averaged curvature and fluctuations of curvature of the mechanical 
manifold associated to the Hamiltonian. This geometrical approach has proved to give 
very accurate estimates in a large number of Hamiltonian physical systems for which 
chaos mainly originates from parametric instability JE^l, that is systems for which the 
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curvature is mainly positive but fluctuating. The method was originally applied to 
derive estimates of the LLE in the thermodynamic limit iV —>■ oo. Later on, this 
approach enabled the analytical computation [Q of the LLE for the HMF model as a 
function of e. In the homogeneous phase (e > e c ) in which the LLE, denoted hereafter 
Ai, vanishes in the limit iV — > oo, it was shown that keeping the leading order in N in 
the ensemble averages of geometric quantities enabled to derive the scaling of Ai as a 
function of N. In this way, the LLE was predicted to scale as iV" 1 / 3 for e > e c . The 
aim of this calculation was to relate the behaviour of the dynamical indicator Ai to the 
occurrence of a second order phase transition at e c in the system. 

In this Letter we wish to apply the same approach to the generalized model ([!]) in 
order to give an analytical prediction for the exponent k. We shall use the derivation of 
the canonical thermodynamics for the system ([]]) presented in reference fl3| . 

Let us first recall the expression for the LLE derived from the geometric approach [§, 
The effective curvature felt by a geodesic is modeled as a Gaussian stochastic 



10, 11 



process whose mean is the average Ricci curvature and the variance its fluctuations. 
Under the ergodic hypothesis, these quantities may be replaced by their microcanonical 
ensemble- averages, denoted respectively kq for the curvature and a 2 , for its fluctuations. 
Then M 



where 



Ai 



A 



A 



4kq 
3A 



(4« /3) 3 + o-y< 



1/3 



and where r is a time scale for the stochastic process estimated as 

T 



2y//Mj + a K o K 



7T 



(3) 



(4) 



(5) 



Using Eisenhart metric, the microcanonical average k q of the mean Ricci curvature kn 
reads k = (k R (9)) where 



k R {9) 



K R {9) 
N - 1 



1 N 

i=l 



d 2 V 
89? 



and 



(6) 



(7) 



As explained in references || |10|, |TTJ we shall assume the equivalence of microcanonical 
and canonical ensembles, leading to identical values of both ensemble-averages of 
observables in the limit A^ — > oo, but to different values of their fluctuations with 



the formula 16 



<* 2 /> „ = + 





\9(f) c ] 


V d(3 J 


9(3 



(8) 
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where (5 = 1/T and kg = 1- This assumption is fully justified by the recent achievement 
that ensemble inequivalence for averages in long-range systems is to be expected only 
close to first order phase transitions [TJ], [|, Therefore, the first step is to express 
the canonical thermodynamics of system (|1|) following reference [|13|]. The constraint 
i ytz j in the potential can be removed for free by defining = b, that is for the 
time being an arbitrary constant. Then the symmetric distance matrix R'^ = r~j a may 
be diagonalized, which enables to use the Hubbard-Stratonovitch transform to evaluate 
the potential part Zy c of the canonical partition function. Using then the saddle-point 
method, one obtains from equations (12) and (20) of reference []13 in the long-range 
case (a < d) and for e > e c (zero magnetization phase) 



InZ 



(3 



N , N 



Vc 



2N^-! " 2 

j=l n=l 



N 



(9) 



The eigenvalues A n of the R' matrix can be easily derived following reference [|19 



Let us first consider the d = 1 case for the sake of simplicity. We remind that 

V* 3 - < io > 

[b if % = J 

where 

Tij = min \i — j + IN\ . (11) 

Therefore R'^ = R'(i — j) = R'(m) where R' is a iV-periodic function. This periodicity 
of the lattice enables to diagonalize R' in Fourier space. Its Fourier transform is 

N 



R'{n) = exp(—i—nm)R'(m) 



N 

m=l 



with the inversion formula 

N 



N ^— ' ^ v iV* 

fe=i 

Then it can be easily shown that R'(i — j) = J2k=i u lk^i u jk where Ujk ■— 
N~ l l 2 exp(-i^rjk) is an element of the unitary matrix of eigenvectors with the following 
expression for the eigenvalues (1 < k < N) 

N 2tt 

X k = J2 eM-i^km)R'(m). (12) 

m=l 

For any d- dimensional lattice, one would get the generalized expression for A n = 
A (ni, . . . , n d ) with 1 < ni, . . . ,n d < N l l d as 

N i/d N i/d f 2n d \ 

An = ■ ■ ■ 2 exp ( nim i ) R '^- ( 13 ) 

mi=l 771^=1 \ *J=1 / 

Coming back to the d — 1 case, we shall take iV even in the following and put iV = 2p. 
As R' is an even function, this implies, for 1 < k < N = 2p, 

X k = b + X k (14) 
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pa 



cos (7rkm/p) 



m=l 



(15) 



\ p is the smallest of the A^'s and is negative. In order to get a fully positive spectrum that 
enables to apply the Hubbard-Stratonovitch transform, we can now shift the spectrum 
by fixing 

p-i 



b := -\ r 



-If 



2 E 



m=l 



(16) 



N is then defined as the maximal eigenvalue 

1 1 

iV = A 2p = 6+ — + 2V — 

' net Z-^/ rn a 



m=l 



(17) 



and it can be easily checked that 

N=2p 

*k = °' 



k=l 



We can now go on with the derivation of the LLE. Using expression @, one gets 



Thus 



That is 



k R {9) 



(k R (0)) c 



<M*)> t 



1 1 N 



cos 



1 



N — 1 



1 N 



N — 1 



1 - 

-E 

AT ^-^ 



1 



d\nZ Vc 

13 3(3 



N 



N ^N-(3\ n 



N-l 

Let us now consider the fluctuations of the curvature. Their canonical average is 



4 d 2 In Z 



N 



Vc 



A, 



(19) 



(20) 



(21) 



(22) 



iV-i d(3 2 N-l^\N-p\ r , 

while the corrective term @, needed to get the fluctuations in the microcanonical 
ensemble, is 



(d{e) c \ l 


\9(k R ) c ] 


V df3 J 


9(3 



(23) 



Actually, the energy density above the critical point is intensive and equal to 1/(2/5) 
(up to a constant term). And one obtains 

1/2 



1 - (5 2 k R ) c (5 2 k R ) 



(24) 
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We shall focus on the high temperature regime and derive the scaling of Ai with 
iV under both the limits of large iV and small /3. We will show later the validity of 
this scaling of Ai with iV in a wider range of /3. The results will be compared with 
the numerical results of references [|14], [HJ for which (3 = 1/9. We shall also compute a 
numerical estimate of the scaling of Ai for (3 = 1/9 using directly equations ([§-§) with 
the full expressions (0-^-0). 

Using ([21]), one gets 

(k R (0)) c ^o +W ^t X2 n (25) 

where En=i K = E*=i ~K + Nb 2 . Using ©, one obtains En=i^ = Np~ 2a + 
2NJ2 P m =i m ~ 2a - If 2a ^ !' the sum of the A n' s is of the order A^ 2 " 20 - Therefore, 
using (H) and (|17D, the scaling of kq with N depends on the value of a as 

{(3N- 1 for < a < 1/2 

p N -2+2a for x/2 < a < i (26) 
/5(lniV)" 2 fora = l. 

Let us now estimate the order of the fluctuations. From fl22|), (p3j) and in the limit of 
vanishing /3, a 2 is equivalent to 2 (N — 1) _1 Af -2 J^ n=1 A 2 , so that 

{A/" -1 for0<a<l/2 
N -2+2a for x/2 < a < i (27) 
(lnA^T 2 fora = l. 

At this stage we can check that the results presented in reference [|[] are effectively 
recovered for the HMF case in the high energy regime where the curvature and 
fluctuations of the curvature were predicted to be of order A/ --1 . Now, considering 
the time scale estimated as (^|), one obtains yfn$ + a K ~ y/b\ that vanishes in the limit 
of large A" while a K j ~ P" 1 ' 2 = 0(1). This last estimate is interpreted in references 



TTJ as the relevant timescale when the fluctuations are of the same order as the 
curvature and does not require the positivity of the curvature, which is effectively the 
case here. It comes out then that r is of order one. Putting this together with estimates 
(PB|) and (^7|) in equations (||) and (|3|) leads to the expression of the LLE, in the limit 
of large N and at the leading order in N, as 

A^ 1/3 for0<a<l/2 

A x (AT, p < 1) oc^ + P 1/6 x <J N~ 2/3+2a/3 for 1/2 < a < 1 (28) 

(lnA0" 2/3 fora = l. 

where one can also check the integrability of the model in the large temperature ft — > + 
limit for all A" [p0[| . If one develops the curvature at higher orders in /3, one obtains 

a power series of the type /? Yln=i (^n/Nj + P 2 Eli \^n/Nj + • • • where it can be 
checked that all the sums of the successive powers of \ n / N are equivalent in the large N 
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limit. Similar expansions are obtained for the fluctuations of curvature. In the large N 
limit, one can thus conclude that equation (23) is valid for small (3, and not only in the 
(3 — > + limit. Formula is hence still valid after replacing /3 1//6 with a function of (5 
whose power series development in (5 is in principle calculable. Moreover, coming back 
to the general ci-dimensional case, the previous procedure can be followed replacing (p"2j) 
by (fT3"|). The sum involves then 

AT TV 1 /'' TV 1 / 1 * 

^~iV^...^W + ... +m r~ NN 1 ~ 2a / d 

n=l mi=l m^=l 

if 2a/d 1. Therefore we can replace a by a/d in all the previous scalings. This 
analysis predicts then the universal exponent Hi cLS 

( 1/3 if < a/d < 1/2 

K = K(a/d) = < 9 (29) 
I Ul-a/d) if 1/2 < a/d < 1. 

Hence k is a function of a/d which is equal to 1/3 for a/d = 0, consistently with 
numerical and analytical results obtained for the HMF model || [7|, §]. It vanishes in 
the limit a/d — > 1 consistently with the intensivity of the LLE for short-range (a > d) 



potentials predicted from thermodynamics |13| . Moreover we predict that for a = d the 
LLE should scale as (In N)~ 2 ^ 3 in the limit of large N and energy density. This agrees 
with the intuitive statement made in reference [14] that the LLE should scale as some 
power of 1/ In N for a — d. 

We shall now compare these estimates with the numerical results presented in 
reference |I[ for d — 1 and in references JR], [15] for d = 1, 2 and 3. In these papers, the 
LLE has been computed numerically for values of N ranging from 5 to 1000 for d — 1 



and from 36 to 3969 for d = 2,3. In references [0, |Tj, [151 , the curves giving the LLE as 
a function of iV in log-log scale were fitted by the functional form aN~ K + bN~ K ~ c with 
c = 1. We do not agree that this form is the correct second order approximation for 
any value of a/d between and 1, nor that c should be equal to 1. Moreover, due to 
obvious numerical constraints, the values of N used are not sufficiently large to discard 
additional finite-iV effects. For the values used in references 0, pl| for d = 2 and 
d = 3 we instead performed a linear fit retaining only the highest values of N, with 
iV between 500 and 4000, as shown in figures || and [| and got the numerical values of 
k plotted in figure |3|. On the other hand, using (04j-iD together with ( pl| - |2l^|2"4"D , we 
can compute the LLE for different values of a and iV and obtain the theoretical curves 
associated to the numerical ones in figures p] and 0. In log-log scale, we have derived 
the best linear fits from the theoretical curves giving the LLE as a function of iV for 
different values of a. It should be mentioned that the results so far obtained for k 
benefit from a better statistics for the largest values of iV compared to the numerical 
results of references O, [T5J, which explains that the theoretical curve be slightly above 



numerical points in figure [| Nevertheles we observe in this figure a good agreement 
between the theoretical and numerical finite- N derivation of k as a function of a/d. 
The agreement is worse for d = 3, compared with the d = 2 case; this may be related 
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to more important finite N effects as d increases for a given value of N . This supports 
the validity of our theoretical predictions. Consequently, we suggest that the so-called 
'universal curves' published in references |[| and [Tj| are plagued by finite- N effects. 



In the asymptotic limit in N, we claim that the correct universal scaling of the LLE 
with N is given by (p9|). This means that there is a sharp change in the dynamical 
behaviour of the model for a/d = 1/2. For < a/d < 1/2, the suppression of chaos 
scales like in the HMF model, which provides then a universal law in this range, with a 
universal exponent k = 1/3. 
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Figure 1. Loglog plot of the LLE as a function of N for different values of a with 
d = 2. The data have been communicated by the authors of reference |Q. At each 
value of a/d, the best linear fit of the logarithm of the LLE as a function of the 
logarithm of N is plotted for N between 400 and 4000. The values of a/d are, from 
top to bottom, 1, 0.95, 0.8, 0.6, 0.4, 0.2 and 0. 
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Figure 2. Same plot as in figure |l] with d = 3. The values of a/d arc, from top to 
bottom, 0.95, 0.8, 0.6, 0.4, 0.2 and 0. 
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Figure 3. Plot of the exponent k as a function of a/d. The bold line corresponds to 
the analytical prediction obtained in the infinite N limit. The thinner line gives the 
exponent obtained by fitting the logarithm of the analytical prediction of the LLE as 
a function of the logarithm of N for N between 500 and 4000. This is to be compared 
with the values of k as a function of a/d deduced from the fits shown in figures ^ (with 
filled squares) and || (with empty circles) . 



